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MEASUREMENTS AND POWER SPECTRA OF 
RUNWAY ROUGHNESS AT AIRPORTS IN COUNTRIES OF 
THE NORTH ATLANTIC TREATY ORGANIZATION 
By Wilbur E. Thompson. 


SUMMARY 


Measurements of runway roughness obtained by a profile- survey method 
are presented. A variety of runway lengths or sections and degrees of 
roughness were measured. The results are presented as elevation pro- 
files of the runways surveyed and in the form of power spectra. 


INTRODUCTION 


Because of concern regarding taxiing loads due to runway roughness, 
various programs have been -undertaken to accumulate data on runway pro- 
files. An initial step in this endeavor was the roughness measurement 
of two runways at Langley Air Force Base, Va. These data are presented 
in reference 1. An additional aid in this study was the roughness meas- 
urement of runways at four commercial airports which have large amounts 
of traffic. These data are presented in the form of elevation profiles 
in reference 2. 

In order to increase the amount of data available for runway rough- 
ness studies, the National Advisory Committee for Aeronautics has agreed 
to process runway-profile data obtained through the Advisory Group for 
Aeronautical Research and Development on runways located in member coun- 
tries of the North Atlantic Treaty Organization. These data are com- 
piled by the surveying method described in reference 1. As a first step 
in this program, measurements have been obtained for 28 runways. This 
report presents measurements of these 28 runways and 6 additional run- 
ways as elevation profiles and power spectral densities. 


SYMBOLS 


number of uniformly spaced points over frequency range at which 
power estimates are desired 


m 
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n 
x 

Ax 
y(x) 

A 

a root-mean-square value of y 


4>(fl) power- spectral-density function, 

0 reduced frequency, 2 tc/'K, radians/ft 


SURVEY OF RUNWAY ELEVATIONS 


The elevation-profile graphs and power-spectral-density graphs for 
the 3k runways evaluated herein are presented in figures 1 to 56 . Run- 
ways 5 R and 5L (figs. 8 and 9) constitute measurements along both the 
right and left sides of one runway. Runways 10-1 and 10-11 (figs. 17 
and 18 ) constitute measurements of two sections or portions of one run- 
way; and likewise, runways 22-1 and 22-11 (figs. 35 and 36 ) represent 
two sections of one runway. On runways 33 and 3 k (figs. 53 and 3k), 
measurements were made each month for a period of nine consecutive 
months . 

The elevation-profile measurements were obtained by the profile- 
survey method described in references 1 and 2 , that is, by means of a 
surveyor’s level, rod, and tape, and by using an interval of 2 feet with 
the exception of the measurements on runway 8 (fig. 1 ^) where an inter- 
val of 0.6 meter (approximately 1.97 feet) was used. As explained in 
reference 1 , the selected measurement interval of 2 feet is one-half 
the shortest wavelength considered and adequately specifies the dis- 
turbance. This 4- foot wavelength corresponds to approximately 35 cycles 
per second at 100 miles per hour. 

The runway elevations for runways 7 and 8 (figs. 13 and 1^) were 
tabulated in meters in the original data and were converted to feet for 
use in this paper. 



number of equally spaced elevations taken over runway 
distance, ft 
space interval, ft 

random function of distance (runway elevation) 
wavelength, ft 
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The measurements were taken along a line parallel to the runway- 
center line, approximately along the path of a main landing gear of an 
aircraft landing or taking off. 

The elevation references are arbitrary; however, the zero reference 
for measurement of several of the runways was close to sea level. 


POWER SPECTRA OF RUNWAY ELEVATIONS 


The power spectral densities which provide a description of the 
frequency characteristics of the runway roughness were computed by the 
method described in reference 3 - Specifically, the computation of the 
power spectra involved the following steps: 

(l) The measured profile data were first prewhitened by the 
relation 


y'(x) = y(x) - 0y(x - Ax) (0 = l) 

This linear operation corresponds to the multiplication in the frequency 
plane by the frequency response function F(fl) given by 


F(ft) = 1 - e~ inAx 


(2) The mean lagged products or autocovariance functions 
determined from 


Rp were 


n-p 


= irH X y,( i y,cl+p (p = ° A ' 2 ’ * • • m) 


q=l 


( 3 ) Initial or "raw" estimates of power were determined by 


Lh 


m 

= ^ e p cos R p (h = 0,1,2, . . . m) 


p=0 
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where 



(0 < p < m) 

(p = 0 or m) 


(4) Final or "smoothed" estimates of power were determined as 

follows : 


= 2 ^0 + 2 Ll 

(h = 0 ) 

*h “ J ^-1 + 2 ^ + t ^+1 

(l S h £ m~- 1 ) 

®h = 2 ^Di-1 + 2 

(h = m) 


( 5 ) The results of the previous step are then divided by 

|f(A)| 2 = 2(1 - cos ft Ax) 

to compensate for the prewhitening operation. 

Evaluation was made with m = 40, which is the number of uniformly 
spaced points over the frequency range investigated. Power- spectral- 
density values are plotted as a function of the reduced frequency 
Plf 

ft = Each point on the power- spectral-density graphs represents the 

A 

average power in a frequency interval ±~ about the value plotted. 

oO 

The cut-off frequency at the high-frequency end of the power curves is 
Radians per foot, with the exception of runway 8 (fig. 15) where the 

cut-off frequency was —2 — radians per foot. 

1.97 

The root-mean- square values were derived from the equation 

* - I 4(0) r m 

r=l 
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since the mean square of runway elevation is represented by the area 
■under the power estimate curves. However, a lower limit of zero fre- 
quency is neither attainable nor desired; therefore, the root -mean- square 
values presented herein are approximations for comparison of average 
roughness within the frequency range investigated. The root mean square 
of runway elevation, c, which is a convenient measure of average rough- 
ness, is presented for each runway on the respective power- spectral- 
density graphs. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., April 28, 1958. 
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Figure 2.- Profile of runway 2. 
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Figure 3 .- Power- spectral-density functions for runways 1 and 2. 
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Figure 7 .- Power- spectral-density functions for runway k 
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Figure 10.- Power-spectral-density functions for runway 5 
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Figure 17*- Profile of runway 10- I. 
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Figure 21.- 

Profile of 

runway 12. 
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Figure 23-- Power- spectral-density functionB for runways 12 and 13. 
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Figure 45-- Profile of runway 26. 
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Figure 1+6 . - Power- spectral-density functions for runways 27 and 28. 
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Figure ^9.- Power- spectral-density functions for runways 29 and 30* 
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Figure 50.- Profile of runway 31 


NACA OT 4303 




NACA TN 4-303 


57 





.01 .1 1 10 

Reduced frequency, q, radians/ft 


Figure 52.- Power-spectral-density functions for runways 31 and 32. 
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(d) March 1957- 
Figure 53*- Continued. 
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(f) May 1957. 


Figure 53 « ~ Continued. , 
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Figure 5^. - Profiles of runway Jk. 
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(fc) January 1957- 
Figure 5^.- Continued. 
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(c) February 1957* 
Figure 54.- Continued. 
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(e) April 1957- 

Figure 5^.- Continued. 
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Figure 55.- Power- spectral-density functions for runway 33* 
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.Figure 55 ,- Continued 
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Figure 55*- Continued 
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Figure 55 .- Concluded 
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Figure 5 6 .- Power- spectral-density functions for runway 3^, 
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Figure 56 .- Continued 
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Figure 56 .- Continued. 
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Figure 56.- Continued 


li ■ liii.i-L, ||, , .. ii j.., , lUfciJlI.Mi'U. I . »«rlfcl li. .. . j. 1 illii 1,. Lj.. I liilLij-HWil: 





